ELLIPTIC CURVES 
AND CONTINUED FRACTIONS 
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Abstract. We detail the continued fraction expansion of the square root of 
the general monic quartic polynomial, noting that each line of the expansion 
corresponds to addition of the divisor at infinity. We analyse the data yielded 
by the general expansion. In that way we obtain 'elliptic sequences' satisfying 
Somos relations. I mention several new results on such sequences. The paper 
includes a detailed 'reminder exposition' on continued fractions of quadratic 
irrationals in function fields. 



1. Introduction 

A delightful 'essay' |l3 by Don Zagier explains why the sequence {Bh)hei , defined 
by _B_2 = 1 , B-i = 1 , Bo = 1, J5i = 1 , B2 = 1 and the recursion 

(1) Bh^2Bh+3 = 5/1+2-8^-1 + Bh+iBh , 

consists only of integers. Zagier comments that 'the proof comes from the theory 
of elliptic curves, and can be expressed either in terms of the denominators of the 
co-ordinates of the multiples of a particular point on a particular elliptic curve, or 
in terms of special values of certain Jacobi theta functions.' 

In the present note I study the continued fraction expansion of the square root of 
a quartic polynomial, inter alia obtaining sequences generated by recursions such 
as Here, however, it is clear that I have also constructed the co-ordinates of the 
shifted multiples of a point on an elliptic curve and it is it fairly plain how to relate 
the surprising integer sequences and the elliptic curves from which they arise. 

A brief reminder exposition on continued fractions in quadratic function fields 
appears as ^ starting at page 1801 below. 

It turns out that many of my results related to sequences a la (QJ are contained 
in the recent thesis ilSj of Christine Swart; for extended comment see ^S] Where 
appropriate, I have added remarks throughout the paper. Michael Somos, see 
had inter alia asked for the inner meaning of the behaviour of the sequences (B;,) , 
above, and of {Ch) defined by Ch-2Ch+2 = Ct-iCh+i+Cl and C_2 = 1, C_i = 1, 
Co — 1, Ci = 1: in the terminology of ^S], of the sequences Somos (5) and 
Somos (4) . More generally, of course, one may both vary the initial values and 
coefficients and generalise the 'gap' to 2m or 2m -|- 1 by studying Somos 2m, 
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respectively Somos 2m + 1 , namely sequences satisfying the respective recursions 

rn m 
i=l i=l 

I show in passing, a footnote on page 23 that a Somos 4 is always a Somos 6, 
while Theorem 121 points out it is always a Somos 5. After seeing jT2], I added a 
somewhat painful proof, Theorem |31 on page[7Hl that it also always is a Somos 8. 
For example, Somos (4) satisfies all of 

Ch-sCh+s, = Ch-iCh+i + 5C^ , 

Ch-iCh+z = —Ch-iCh+2 + 5ChCh+i , 

Ch-iCh+i — 25C;i_iC;i+i — 4C,^ . 

In the light of such results one can be confident, see '>\i.'6\ page [73 that in general 
if {Ah) is a Somos 4, equivalently an elliptic sequence satisfying Ah-2Ah+2 = 
WiAh-iAh+i - W1W3AI, then for all m both 

WiW2Ah-„iAh+m+l = WmWm+lAh-lAh+2 - Wm-lWm+2AhAh+l 

and 

W^Ah-^Ah+„i = W^Ah-lAh+l - Wm-lWrn+lAl . 

2. Continued Fraction Expansion of the Square Root of a Quartic 

We suppose the base field F is not of characteristic 2 because that case requires 
nontrivial changes throughout the exposition and not of characteristic 3 because 
that requires some trivial changes to parts of the exposition. We study the continued 
fraction expansion of a quartic polynomial D G ¥[X] ; where D is not a square. 
Set 

(2) C :Y'^ ^ D{X) -.^ {X^ + ff + Av{X -w), 

and for brevity write A = X^ + f and R ~ v{X — w) . For /i = , 1 , 2 , ... we 
denote the complete quotients of Yq by 

(3) Yh = {Y + A + 2eh)/vh{X-Wh), 

noting that the Yh all are reduced, namely degY/i > but degF^ < 0. The 
upshot is that the continued fraction expansion of Iq has typical line, line h : 

Y + A + 2eh ^ 2{X + Wh) _ Y + A + 2eh+i 
Vh{X - Wh) Vh Vh{X - Wh) 

Thus evident recursion formulas, see at page|Sni yield 

(4) f + eh + e/i+i = -wl 

and -VhVh+i{X - Wh){X -Wh+i) ^ (Y + A + 2eh+i){Y + A + 2eh+i). Hence 

(5) VhVh+i{X - Wh){X - Wh+i) = -4(^2 + / + eh+i)eh+i + 4i; (X - w). 
Equating coefficients in l^, and then dividing by —Ach+i-, we get 

X'^- -4:eh+i = Vhi'h+i; 

X^- v/eh+i^Wh+Wh+i; 

f + eh+i+vw/eh+i=WhWh+i. 
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The five displayed equations immediately above readily lead by several routes to 

(6) ehCh+i = v{w - Wh) ■ 

For example, apply the remainder theorem to the right hand side of ((HJ after noting 
it is divisible by X — Wh, and recall I^J . 

Proposition 1 (Adams and Razar Denote the two points at infinity on the 
elliptic curve ^ by S and O , with O the zero of its group law. The points 
Mh+i '■= {wh, efi — e/i+i) all lie on C . Set Mi = M , and Mh+i —: M + Sh- Then 
Sh = hS . 

Proof. The points Mh+i lie on the curve C : Y"^ ~ D{X) because 

{eh ~ eh+if - {wl + ff = {eh - {en+i + wl + /)) {{eh + wl + f)- ch+i) 

= -AchCh+i = ■iv{wh - w) . 

The birational transformations 

(7) X={V-v)/U, Y = 2U-{X^ + f); 
conversely, 

(8) 2U = Y + X^ + f, 2V = XY + X'^ + fx + 2v, 

move the point S to (0,0), leave O at infinity, and change the quartic model to a 
Weierstrass model 

(9) W -.V^ -vV ^U^ - fU^ + vwU . 

Specifically, one sees that U{Mh+i) = —Ch+i, and V{Mh+i) — v — Whet+i- We 
also note that U{-Mh+i) = -e/i+i , V{-Mh+i) = WhCh+i- 

To check S + {M + Sh-i) = M + Sh on W it suffices for us to show that the 
three points (0,0), {—eh,v — Wh-iCh) , a,Tid {—eh+i,Wheh+i) lie on a straight line. 
But that is {v — w^-ieh) / eh — Wh ■ So w^-i + Wh— v/^h proves the claim. | 

3. Elliptic sequences 

Theorem 2. Let {Ah) he the sequence defined by the 'initial' values , Ai and 
the recursive definition 

(10) Ah-iAh+i = ChAl . 

Then, given Aq , Ai , A2 , A3, satisfying HfO|l . the recursive definition 

(11) Ah-2Ah+2 = y-'Ah-iAh+i + i;2(/ + ^2)^2 
defines the same sequence as does (|f 0|) . Just so, also 

{12) Ah-2Ah+3 = -v^{f + w^)Ah-iAh+2 + v'^{v + 2w{f + w^))AhAh+i 
defines that sequence. 
Proof. By (0 we obtain 

Ch^ielch+i = v'^{w - Wh-i){w - Wh) 
= v'^{wh-iWh-w{wh-i + Wh) + 'w^) = ((/ + e/i + vw/ch) - w ■ {v/ch) + w^) . 
Thus 

(13) eh-ie\eh+i = {ch + {f + w"^)) . 
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However, Ah-iAh+i = ShAf entails 

and so Ah-2Ah+2 = eh-iCheh+iAh-iAh+i, which is 

(14) Ah-2Ah+2 = eh~ieleh+iAl . 

On multiplying ((T^ by A^ we obtain ifTTl) . 

Similarly ((TUJ) yields Aft+iyl,iAh+2 = ehCh+iAlAl^^, and so 

(15) Ah-iAh+2 = Gheh+iAhAh+i ■ 
It follows readily that 

(16) ^/!-2^/i+3 = e/i_iefje^^;^e,i+2^h • 
Moreover, implies that 

eh-ielel^ieh+2 = i''' (ehe/^+i + (/ + w^)(e/i + e/,+i) + (/ + w^)^) . 

However, by Q we know that v'^ (^Ch + ei^+i + f + w'^) ~ v^{w^ — wf^) . Here 
v{w - Wh) = e^e/i+i and v{w + Wh) = -v{w ~ Wh) + 2vw — -ehCh+i + 2vw. So 

(17) eh-ielel^j^eh+2 = (-(/ + 'u?)eheh+i + + 2vw{f + w^)) , 

which immediately allows us to see that also (|12|l yields the sequence {Ah) ■ | 

The extraordinary feature of the identities H13(l and H17|l is their independence 
of the translation M : thus of the initial data wo , ifo , and eo . 

3.1. Two-sided infinite sequences. It is plain that the various definitions of 
the sequence {Ah) encourage one to think of it as bidirectional infinite. Indeed, 
albeit that one does feel a need to start a continued fraction expansion — so one 
conventionally begins it at Iq i one is not stopped from thinking of the tableau 
listing the lines of the expansion as being two-sided infinite; note the remark at the 
end of S 16.11 page l8UI In summary: we may and should view the various sequences 
{eh) defined above, as two-sided infinite sequences. 

3.2. Vanishing. If say Vk ^ 0, then line k of the continued fraction expansion of 
Yq makes no sense both because the denominator Qk{X) := Vk{X — Wk) of the 
complete quotient Yf. seems to vanish identically and because the alleged partial 
quotient := 2{X + Wfe)/wfe blows up. 

The second difficulty is real. The vanishing of entails a partial quotient 
blowing up to higher degree. We deal with vanishing by refusing to look at it. 
We move the point of impact of the issue by dismissing most of the data we have 
obtained, including the continued fraction tableau, and keep only a part of the 
sequence (e^). That makes the first difficulty moot.* 

Remark. There is no loss of generality in taking fc = 0. Then, up to an irrelevant 
normalisation, Yq = Y + A. If more than one of the Vh vanish then it is a simple 
exercise to confirm that the continued fraction expansion of Yq necessarily is purely 
periodic, see the discussion at page|S21 If does not have a periodic continued 

*In any case, the first apparent difiiculty is just an artifact of our notation. If, from the start, 
we had written Qh = v^X + yji , as we might well have done at the cost of nasty fractions in our 
formulas, we would not have entertained the thought that = entails = 0. Plainly, we 
must allow Wfc = yet visW}^ ^ . 
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fraction expansion then there is some ho, namely /iq = 0, so that, for all h> ho, 
line h of the expansion of Yq does make sense. 

Except of course when dealing explicitly with periodicity, we suppose in the 
sequel that if = then fc = ; we refer to this case as the singular case. 

3.3. The singular case. We remark that in the singular case the sequence {eh)h>i 
defines antisymmetric double-sided sequences (Wh) , that is with W-h = —Wh , by 
Wfi-iWh+i = ehW^ and so that, for all integers h, m, and n, 

m) Wh^mWh+^mW^ + Wn-hW„+hWl + Wm-nWrn+nWi = . 

Actually, one may find it preferable to forego an insistence on antisymmetry in 
favour of rewriting I^T)l less elegantly as 

(18) Wh-mWh+mW^ = Wh-nWh+nWl - Wm-nWm+nW^ , 

just for h > m > n. In any case, (|18|) seems more dramatic than it is. An easy 
exercise confirms that, if Wi = 1 , (|18|l is equivalent to just 

(19) Wh^mWh+m = WlWh^iWn+l - W^.iW,n+lWl 

for all integers h > m. Indeed, 119(1 is just a special case of ((18|l . However, 
given 1)19(1 . obvious substitutions in ((18(1 quickly show one may return from 1(19(1 to 
the apparently more general ((18(1 . 

But there is a drama here. As already remarked in a near identical situation, 
the recurrence relation Wh-2Wh+2 — W2Wh-iWh+i — WiW^W^ , and five or so 
initial values, already suffices to produce (Wh) ■ Thus ((19(1 for all m is apparently 
entailed by its special case m = 2. 

I can show this directly^, by way of new relations on the eh, for m = 3. But 
the case m = 4 already did not seem worth the effort. Whatever, my approach 
gave me no hint as to how to concoct an inductive argument leading to general m . 
Plan B, to look it up, fared little better. In her thesis ^2], Rachel Shipsey shyly 
refers the reader back to Morgan Ward's opus ^HI; but Ward does not comment 
on the matter at all, having defined his sequences by ((19(1 . Well, perhaps Ward 
does comment. The issue is whether ((19(1 is coherent: do different m yield the 
one sequence? Ward notes that if cr is the Weierstrafi cr -function then a sequence 
{a{hu) / a{u)'^ ) satisfies (119(1 for all m. Whatever, a much more direct argument 
would be much more satisfying. 

Proposition H shows that certainly Wh-2Wh+2 = WiWh-iWh+i - WiWsW^ 
for h = 1, 2, . . . , in which case ((1911 apparently follows by arguments in |E| and 
anti-symmetry; ((18(1 is then just an easy exercise. 

The singular case is initiated by wi = 4w, wi = w, ei = 0, 62 = — (/ + w^). For 
temporary convenience set x — v/{f + w'^) . From the original continued fraction 

tpiainly eft_2e2_^e3 e2_|_^e^+2 ■ = ^'^(eh-i + (/ + «'^)) (e/i+l + (/ + w^))eh- Now notice 
that {eh-ieh + eheh+i)eh = v{w -w^-i+w -Wh)eh = 2vweh-v^ and recall that e^-ieleh+i = 
(fift "I- (/ "I- w^)) ■ The upshot is a miraculous cancellation yielding 

eh-2eLl4eh+ieh+2 ■eh=v'^{{f + w'^fel + v(v + 2w{f + w^))eh) 

and allowing us to divide by the auxiliary . Thus the bottom line is 

Ah-3Ah+3 = v*{if + w^Ah-iAh+i + v(v + 2w(f + w^))AI) , 

which is Ah_3Ah+3 = W^A^-iA^+i - W2W4AI . 

*For additional remarks, and a dissatisfying proof for the case m = 4, see il5.2l at page I78l 
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expansion of Y + A or, better, the recursion formulas of page 1701 we fairly readily 
obtain V2 = l/x, W2 — w~x, 63 = ~x{x+2w) , 64 — v{x'^(x+2w) — v)/x'^{x+2wY ■ 

We are now free to choose, say Wi — 1, W2 = v , leading to W3 — —v'^{f + w^), 
Wi = -v'^{v + 2w{f + w'^)) , = + 2w{f + w^)) - (/ + ^2)3) , . . . . 

That allows us to notice that (|12|) apparently is 

and that Hll() of course is 

Ah-2Ah+2 = WiAh-iAh+1 - W1W3AI . 

Given that also A/i_3yl/i_|_3 = W^Ah-lAfl-^-l — W2WiA\ , it is of course tempting to 
guess that more is true. Certainly, more is true in the special case {Ah) — {Wt) ^ 
that's the point of the discussion above. Moreover, the same 'more' is true, see for 
example ^| Theorem 8.1.2, p. 191], for sequence translates: thus [Ah] = [Wh+k) ■ 

3.4. Elliptic divisibility sequences. Recall that in the singular case and for 
= 1 , 2 , ... the — e/i are in fact the U co-ordinates of the multiples hS of the 

point S = (0, 0) on the curve V"^ -vV = U^ ~ fU"^ + vwU . 

Suppose we are working in the ring Z = Z[/, v, vw] of 'integers'. If gcd(w, vw) — 
1 , so the exact denominator of the 'rational' w is v, then our choices Wi — 1 , 
14^2 = V lead the definition Wh-iWh+i — enW^ to be such that is always the 
exact denominator of the 'rational' Ch ■ It is this that is shown in detail by Rachel 
Shipsey jj^l- In particular it follows that {Wh) is an elliptic divisibility sequence as 
described by Ward |lf)j . A convenient recent introductory reference is Chapter 10 
of the book [2]. 

Set hS = {Uh/W^,Vh/W^), thus defining sequences (Uh) , (K) , and (Wh) 
of integers, with Wh chosen minimally. Shipsey notices, provided that indeed 
gcd{v,vw) = 1, that gcd{Uh,Vh) = Wh-i and Wh-iWh+i = -Un- Here, I have 
used this last fact to define the sequence {Wh) ■ 

Starting, in effect, from the definition (|18|l . Ward shows that with Wq — 0, 
Wi — 1, and VF2IVF4, the sequence {Wh) is a divisibility sequence; that is, if 
a\b then Wa\Wb- A little more is true. If also gcd(W3,W4) — 1 then in fact 
gcd(Wa, Wf,) = Wgcd(a,b) ■ the other hand, a prime dividing both W3 and W^ 
divides Wh for all /i > 3. 

3.5. Periodicity. Suppose now that the sequence {Wh) is periodic. From the 
continued fraction expansion and, say, 0, we find that v = Q (but w' ~ vw 7^ if 
our curve is to be elliptic) is the case of the continued fraction having quasi-period 
r = 1 and the divisor at infinity on the curve having torsion m — 2. Just so, 
f + w-^ = , thus T4^3 = , signals r = 2 and m = 3 , and x -I- 2?« = , or W4 = , is 
r = 3 and m — -i. And so on; for more see 8; . In summary, to > is minimal with 
Wm = if and only if the continued fraction expansion of, say, Z — ^{Y + A) has a 
minimal quasi-period of length r = to — 1 . Note that by observations summarised 
at ^GA[ page 1821 below, the minimal period of the continued fraction expansion of 
Z then is m — 1 or 2to — 2 if to is even, and m — 1 if to is odd. This holds ^ for 



^ This fact is not well understood. For example, over the finite field ¥p of p elements, p an odd 
prime, it is plain and well known that a quasi-period of length r entails there is a period of length 
rt for some divisor t of p — 1 . Indeed, if we are expanding a general quadratic irrational function, 
then this is the best one can assert. If, however, the function has integral — thus polynomial — 
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any base field. Over Q, recall p^, the only possibilities for m are 2, 3, . . . , 10, or 
12; by happenstance, the cases period 9 or 11 do not occur (see ^Hl or [H]). 

However, Ward |16| shows and one fairly readily confirms that precisely the 
periods 1, 2, 3, 4, 5, 6, 8, or 10 are possible for an integral elliptic divisibility 
sequence defined by (|19|l . A convenient summary and explanation is given by 
Chistine Swart PI §4.6]. 

Aside: It has been suggested in my hearing that 'Mathematics is the study of degeneracy', so 
the following warrants careful consideration. In the singular case we have ei = and then 
the recursion e(,_ie^e(,^i = v^i^e/^ + {/ + ^^)) and ^2 = — (/ + w^) forces eg ■ 0^ = —v^ in 
the case h = \. In a context in which d = and w' = vw ^ passes without comment 
this is no great matter. However, we must also cope with W—\W\ = coWq , and with both 
W—2W0 = e—iW^i and W—sW— 1 = e_2W^2 ^^'^ now is more difficult to reconcile the two- 
sided continued fraction expansion and the two-sided anti-symmetric singular elliptic sequence. 
The periodic case seems particularly problematic. Noticeably, Christine Swart declares her 
elliptic sequences as undefined beyond a ; I've chosen to be vague. 

4. Examples 

4.1. Consider the curve C : — {X^ — 29)^ — 4 • 48{X + 5) ; here a corresponding 
cubic model is £ : + A8V ^ + 29U^ + 2A0U . Set A = X^ -29. The 
first several preceding and succeeding steps in the continued fraction expansion of 

X-2 Y + A + 32 
8/3 " 16(X-f 2)/3 
X -1 Y + A + 24 
~~6 12(X + 1) 

X-3 Y + A+16 
~ 4(X + 3) 

X-3 Y + A + 2A 
~ 8{X + 3) 

X-1 Y + A + 32 
~ 6{X + l) 

X-2 F + yl + 18 
16/3 ~ 32{X + 2)/3 



where elegance has suggested we write 'line h' as short for 'line —h'. 

The feature motivating this example is the six integral points (—2, ±7) , (—1, ±4) , 
and (—3, ±4) on C. With Mc — (—3,4) and Sc the 'other' point at infinity these 
are in fact the six points Mq + hSc for h = —3, —2, —1, 0, 1, and 2. 

Correspondingly, on £ we have the integral points M + 2S — (—16,-32) and 
M -2S ^ (-16,-16), M - S = (-12,-36) and M + S = (-12,-12); here 
M = Mg ^ (-8, -24) ; S* = S'f = (0, 0) . Of course £ is not minimal; nor, for that 



Yo = {Y + A+ 16)/8(X -I- 3) are^ 
line 3 : 



line 2 : 
line 1 : 
line : 
line 1 : 
line 2 : 



Y + A- 


-f 18 


16(X + 


2)/3 


Y + A 


+ 32 


12{X 


+ 1) 


Y + A 


-f 24 


4(XH 


K3) 


Y + A 


+ 16 


8(XH 


K3) 


Y + A 


+ 24 


6(XH 


HI) 


Y + A- 


■f 32 


?,2{X + 


2)/3 


y + aa 


- 18 



9(3X + 10)/8 



trace, a fortiori if it is a quadratic irrational integer element, then if its quasi-period is of even 
length it must already be a period, and if it of odd length r , then at any rate 2r is a period length. 
^Here my choice of ijq = 8 is arbitrary but not at random. 
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matter was C . In fact the replacements X ^ 2X + 1, Y ^ AY yield 
(20) Y"^ = [X"^ + X -if -A-Q{X + i) , 

correctly suggesting we need a more general treatment than that presented in the 
discussion above. It turns out to be enough for present purposes to replace e/j 4e^ 
obtaining 

• ■ ■ , e_3 = |, e_2 = 4, e_i =3, eo = 2, ei 3, 62 = 4, 63 = |, ... . 
Then = 1 , Ai = \ and 

Ah-iAh+i = ehA\ 

yields the sequence . . . , ^-4 = 253^ A_3 = 2^3^, A^2 = 2^3, 2, ^0 = 1, 

Ai = 1, A2 = 3, ^3 = 2^32, Ai = 2^3^, .... Notice that we're hit for sixH 
by increasingly high powers of primes dividing 6 appearing as factors of the Ah . 
However, we know that p2l) derives from (|17|l . With the original e^s divided by 4 
that yields 

Ah^2Ah+?, = 6^Ah-iAh+2 + Q^AhAh+i . 
Remarkably, one may remove the effect of the 6 by renormalising to a sequence 
(Bh) of integers satisfying 

Bh-2Bh+3 = Bh^iBh+2 + BhBh+i ■ 

Specifically, ... , B_4 = 3, S_3 = 2, S_2 = 1, B-i = 1, Bq = I, Bi = 1, 
B2 = I, B3 = 2, Bi = 3, = 5, Be = 11, Bj ^ 37, Bg = 83, and 
the sequence is symmetric about B ~ Q. Interestingly, the choice of each Bh as 
a divisor of Ah is forced, in the present case by the data Aq = Ai = 1 and the 
decision that the coefficient of BhBh+i be 1 . Of course it is straightforward to 
verify that Ah-2Ah+z is always divisible by 6'^ and Ah-iAh+2 always by 6. For 
a different treatment see i)4. 41 below. 



4.2. Take v — ±1 and f + uP' = 1. Thus eh-ie\eh+i = + 1 and so eo = 1, 
ei = 1 yields the sequence 2, 1, 1, 2, 3/4, 14/9, . . . , of values of eh- As 
explained above, with Cq — 1 and Ci = 1, the definition Ch-iCh+i — ^hCh yields 
the symmetric sequence . . . , 2 , 1 , 1 , 1 , 1 , 2 , 3 , 7 , 23 , 59 , . . . , of values of Ch 
satisfying the recursion 

Ch-2Ch+2 = Ch-iCh+i + Ch ■ 

Plainly, one can get four consecutive values 1 in a sequence {Ch) as just defined 
only by having two consecutive values 1 in the corresponding sequence (ch) ■ 

Set Y^ ^ A^ + Av{X - w) , where A = X^ + f . With Z = ^{Y + A) , we have 
ZZ = —v{X — w) and Z + Z — A. Thus the consecutive lines 

iX + b)-^+' 



X-b ' ' X-b 
Z + 1 , Z + 2 



-{X-c) ' ' -{X-c) 
entail 

/ + 2 = -6^ / + 3 = -c^ and b + c^v, be ^ f + vw + 1 ^ 



" My remark is guided by knowing that V'^ + UV + 6V = + 7U'^ + 12U is a minimal model 
for £ , and noticing that gcd(6, 12) = 6. 
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which is ti = ±1, b = ±1, c = 0, / = —3, and w = ±2. Up to X < X, the 

sequence (Ch) is given by the curve C : ^ {X'^ — 3)^ + 4(X — 2) and its points 
Mc + hSc , Mc — (1, 0) , Sc the 'other point' at infinity; equivalently by 

£ ■.V'^ -V = U^ + + 2U with M = (-1, 1), S = (0, 0) . 

Indeed, M + S = (-1, 0) , M + 25* = (-2, 1) , M + 2S = (-3/4, 3/8) , . . . . Note 
that it is impossible to have three consecutive values 1 in the sequence (eh) if also 
V — ±1 , except for trivial periodic cases, so the hoo-ha of the example at i)4.1l above 
is in a sense unavoidable. 

4.3. Remarks. The two examples get a rather woolly treatment at ^31 and its 
preceding discussion, see fSj for context. Before seeing ^3] I had also remarked that 
"the observation that a twist V'^ — vV — dU^ — fU^ + vwU becomes ~ dvV — 

~ fU^ + dvw hy U <~ dU , V ^ dV allows one to presume v = ±1. A 
suitable choice of eo, ei and Aq, Ai should now allow one to duplicate the result 
claimed in ^14 in somewhat less brutal form." Namely, one wishes to obtain elliptic 
curves yielding a sequence (Ah) with nominated A_i , Aq , Ai , A2 and such that 
Ah-2Ah+2 = K,Ah-2Ah-2 + ^Aii-2Ah-2 ; only the cases k not a square are at 
issue. In fact, this issue is dealt with by Christine Swart at [131 p. I53ff] in more 
straightforward fashion than I had in mind. In very brief, if Ah-lAfl-^-l — e^A^, 
then 

Bh = K^''^''+^^Ah entails Bh^iBh+i = kchBI , 
and so Bh-2Bh+2 = {1/ Bh-2Bh-2 + (X/ k^)Bi,.2Bi,-2 ■ 

4.4. Reprise. It seems appropriate to return to the example of i)4.1l so as to 
discover the elliptic curve giving rise to {Bh) = (...,1,1,1,1,1,...), given that 
Bh-2Bh+z = Bii-iBh+2 + BhBh+i- Recall we expect the squares of the integers 
Bh to be the precise denominators of the points M + hS on the minimal Weierstrafi 
model W of the curve; here M is some point on that model and S = (0, 0) . 

Suppose e_2 , e_i , eo , ei , 62 supply the five integer co-ordinates yielding i?_2 , 
B-i , Bq , Bi, B2 ■ Of course no more than two of these can be 1 so we must 
have 

eoB-iBi = cqBq , eiBoB2 = eiBf , ^e2BiB3 — e2-B| , 
since of course the recursion for (-B/J entails B^ — 2. Suppose in general that 

ChBh-iBh+i = ChBf^ . 

Then the identity l(T7|l . 

eh-ielel^j^eh+2 = (-(/ + w'^)eheh+i + v"^ + 2vw{f + w^)) , 
and Bh-2Bh+i = Bh^iBh+2 + BhBh+i entail 

Ch-ic\c\^^Ch+2 = -v^{f + w^yhCh+i ^v^{v + 2w{f + w^)). 
Thus CfiCh+i = kv , say, is independent of h and we have 

fc2 = -(/ + w2) and k^ + 2wk^-v = 0. 
Note that if / + , or 2?ii and v, are integers, also k must be an integer. Also, 



(21) 



6961 = vk and 6162 = 2vk 
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Remark. However, e/i_ie^e/i+i = v^Ch + + uP) implies 
k^Bh-2Bh+2 = ChBh^iBh+i + if + w^)Bl , 

without the coefficients necessarily being independent of /i. In particular, — 
— (/ + w^) entails cg = cq = 2fc^ and ci — ei — 3fc^ . 

On the other hand, the identity © now reports that k = w—wq and 2k = w—wi . 
By ^ we then have 

wq + wi ^ v/ei = 2w — 3fc 
whilst by Q we see that / + cq + ei = — (w — k)'^ , / + ei + e2 = — (w — 2fc)^ , so, 
recalling that 62 — 2eo, also eo = (2w — 3fc)fc. Hence 

(22) {2w - 3k)k = 2k^ and so 2w = 5k . 
In summary, we then can quickly conclude that also 

(23) v = 6k^ , 4/ = -29fc2 , and 2wo = 3/fc . 

The normalisation fc = —2 retrieves the continued fraction expansion at M.ll on 
page 123 As shown at on page |HS1 the corresponding minimal WeierstraB model 
is V'^ + UV + 6V = + + 12U; and M = (-2, -2) is a point of order two. 

5. SoMOS Sequences 

5.1. Christine Swart 's thesis |]Jj. Much of the work reported by me here is 
sui generis with my original object having been little more than to make explicit 
the ideas of Adams and Razar ^ and to rediscover elliptic torsion surfaces by a 
seemingly new method, see |S] and its references. Eventually, I learned of Michael 
Somos's sequences, see and realised how they arise from my data. I had sort of 
heard of Christine Swart's work from Nelson Stephens in 2003 but, fortunately as 
it turned out**, did not have access to her thesis ^Hj until very recently when this 
paper was already essentially complete; see [S]- 

Christine Swart's discussion of the interrelationship between elliptic sequences 
and elliptic curves is more detailed and complete than mine. Among many other 
things, she is careful to recognise that the formulas neither know nor care whether 
the given elliptic curve is in fact elliptic: thus, for example, also my quartics may 
have multiple zeros. If so, extra comment — mostly, quite straightforward — is 
required at a number of points above: but is neglected by me. Further, Christine 
Swart reports inter alia that Nelson Stephens (personal communication to her) 
had noticed identities equivalent to (|13(l and dJ, both at page El these are her 
Theorems 3.5.1 and 7.1.2 [13 p. 29 and p. 153]. 

5.2. Much more satisfying? I can in fact show that an elliptic sequence (Ah) 
also is given by A;i_4A;i+4 = W^Ah^iAh+i — W3W5AI. However, I consider my 
argument as typical of the sort of thing that gives mathematics its bad name: and 
regret to have to admit that this sort of nonsense seemingly does generalise to 
proving that a Somos 4 also always is a three-term Sonios (4 + n) . 

Specifically, we know that 

(24) Ah-2Ah+2^WiAh-iAh+i-WiW3Al, by de&iition; 

(25) Ah-3Ah+3 = W^Ah-iAh+i - W2WiAl , see footnote pageESl 



**It was fun to puzzle out not just the answers to questions, but also to attempt to guess what 
the questions ought to be. 
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and intend to show that 

(26) Ah-Uh+i - W^Ah-iAh+i - W3W5AI . 

Notice that, in particular, W1W5 = W^W2Wi - WiWsWi ; so because M^i = 1, 
Wi = WiWaWi -W3W5. Now observe that 

Ah^Uh+2Ah^2Ah+i - iWiAh-2Ah - W2WUl-i)iWiAhAh+2 - W2WiAl^^) . 

In the product on the right, the first term is 

{WiW^WiAl - W3W5Al)Ah-2Ah+2 

and half of it contributes half of (|26|) . Similarly, half the final term of the product, 
thus of 

W^Ah-iAh+i ■ WiAh-iAh+i = WlAh-iAh+i{Ah-2Ah+2 + W3AI) , 

provides the other half of (|26|) . Thus it's ugly but true that we have proved that 
(EH) holds if and only if W3W4 = or 

WiAh-2AlAh+2 ~ W2W3iAh-2AhAl+, + Al_^AhAh+2) + W^Ah-iAlA^+i = . 

I now compound this brutality by fiercely replacing the two occurrences of Ah-2 
by the evident relation 

Ah-2 = {W^Ah-iAh+i ~ WiW3Al)/Ah+2 ■ 

That necessitates our then multiplying by Ah+2 ■ Fortunately, we can compensate 
for this cruelty by dividing by A^ ■ We are left with needing to show that 

(27) WiAh^iAhAh+iAh+2 - WiW3AlAh+2 - M^lT^a^h-i^^+i 

- M^2l^3^Li^'+2 + W2WiAlAl+^ + W4Ah-iAhAh+iAh+2 = . 

It's now natural to despair, and to start looking for a Plan B. However, one might 
notice, on page [731 that W4 = —v'^iyV + 2w{f + w^)) ; and W2 = v. Moreover 
W3 = —v^{f + w^) . Thus, conveniently. 

Hence, just as our result is trivial if W3W4 = 0, so also it is trivial if W2 — 0. All 
this is a sign that we may not as yet have made an error. We may divide (|27|l by 
W2W3 . Better yet, let's also divide by Af^Af^_^_-^ by using the definitions 

Ah-iAh+i = ChAl , whence also Ah-iAh+2 = ehCh+iAhAh+i . 
Then all that remains is a confirmation that 

(28) 2vweheh+i - v'^{eh + Ch+i) - e\e\^-^ - v^{J + w^) = . 

However ((T^ . page[^ is ehCh+i ~ v{w — wt) , while eh + eh+i = — / — is Q), 
page [701 Astonishingly, the claim H28I) follows immediately. 

Theorem 3. // (Ah) is a Somas 4 then it is a Somas 8 of the shape 

Ah-iAh+i = nAh-iAh+i + XAl . 



Proof. Given the argument above, it suffices to note that Christine Swart [151 
p. 153j(f] shows that any Somos 4 is equivalent to an elliptic sequence. | 
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6. Rappels 

6.1. Continued fraction expansion of a quadratic irrational. Let Y = Y(X) 
be a quadratic irrational integral element of the field ¥{{X^^)) of Laurent scries 

oo 

(29) f-hX^^, some d G Z 

h=-d 

defined over some given base field F; that is, there are polynomials T and D 
defined over F so that 

(30) Y"^ =T{X)Y + D{X). 

Plainly, by translating F by a polynomial if necessary, we may suppose that 
deg £> > 2 deg T + 2 , with deg D = 2.g + 2 , say, and deg T <g; then AegY = g + l. 
Recall here that a Laurent series (|29|) with fd^O has degree d. 

Set Yq = (Y + Po)/Qq where Pq and Qo are polynomials so that Qo divides 
the norm (Y + Po){Y + Pq); notice here that an F[J'('] -module {Q,Y + P) is an 
ideal in ¥[X, Y] if and only if Q\{Y + P){Y + P) . 

Further, suppose that degYb > and degFo < 0; that is, Yq is reduced. Then 
the continued fraction expansion of Yq is given by a sequence of lines, of which the 
h-th is 

(31) Yh:^{Y + Ph)/Qh = ah-iY + Ph+i)/Qh; in brief F,, = a,, - B,, . 

Here the polynomial is a partial quotient, and the next complete quotient Y/i+i is 
the reciprocal of the preceding remainder —{Y + Ph+i)/Qh ■ Plainly the sequences 
of polynomials (Ph) and (Qh) are given by the recursion formulas 

(32) P^, + Ph+i + iY + Y) = ahQh and YY + {Y + Y)Ph+i + = -QhQh+i ■ 

It is easy to see by induction on h that Qh divides the norm {Y + Ph){Y + Ph) ■ 
We observe also that we have a conjugate expansion with /i-th line 

(33) Bh := {Y + Ph+i)/Qh = ah - (Y + Ph)/Qh , that is, Bh = Oh - Yh ■ 

Note that the next line of this expansion is the conjugate of the previous line of its 
conjugate expansion : conjugation reverses a continued fraction tableau. Because 
the conjugate of line is the last line of its tableau we can extend the expansion 
forming the conjugate tableau, leading to lines h — \ , 2 , ... 

{Y + P^h+i)/Q-h=a-h-{Y + P-h)/Q-h\ that is, = a_,, - F_,, . 

Plainly the original continued fraction tableau also is two-sided infinite and our 
thinking of it as 'starting' at Yq is just convention. 

6.2. Continued fractions. One writes Iq = [^0 : j ^2 , • • • ] , where formally 

(34) [ao , ai , 02 , . . . , a;i] = flo + 1/ [ai , 02 , . . . , a/i-i] and [ ] = oo . 
It follows, again by induction on /i, that the definition 

Xh Xh-l 

Uh yii-i 

entails [ oq , ai , 02 , . . . , a/j ] — Xh/uh- This provides a correspondence between 
the convergents Xh/yh and certain products of 2 x 2 matrices (more precisely. 





: 
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/ ah 




U 
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between the sequences (xh) , (yh) of continuants and those matrices). It is a useful 
exercise to notice that Yq = [ao , ai , . . . , ah , Yh+i ] impHes that 

Yh+i = -{uh-iY - Xh-i)/{yhY - Xh) 

and that this immediately gives 

(35) YiY2---Yh+i = {-lf{xh^yhY)-\ 

The quantity — deg{xh — yhY) = deg yh+i is a weighted sum giving a measure 
of the 'distance' traversed by the continued fraction expansion to its {h + l)-st 
complete quotient. Taking norms yields 

(36) {xh - yhY){xh - yhY) = {-if+^Qh+i • 

6.3. Conjugation, symmetry, and periodicity. Each partial quotient ah is the 
polynomial part of its corresponding complete quotient Yh . Note, however, that 

the assertions above arc independent of that conventional selection rule. 

One readily shows that Yq reduced, to wit deg lo > and deg Yq < 0, implies 
that each complete quotient Yh is reduced. Indeed, it also follows that deg Bh > 0, 
while plainly dcgi3;i < since —Bh is a remainder; so the Bh too are reduced. In 
particular ah , the polynomial part of Yh , is also the polynomial part of Bh ■ 

Plainly, at least the first two leading terms of each polynomial Ph must coincide 
with the leading terms of Y — T . It also follows that the polynomials Ph and Qh 
satisfy the bounds 

(37) degPh= g + 1 and deg Qh < 9 ■ 

Thus, if the base field F is finite the box principle entails the continued fraction 
expansion of Yq is periodic. If F is infinite, periodicity is just happenstance. 

Suppose, however, that the function field ¥(X, Y) is exceptional in that Yq , 
say, does have a periodic continued fraction expansion. If the continued fraction 
expansion of Yq is periodic then, by conjugation, also the expansion of Bq is 
periodic. But conjugation reverses the order of the lines comprising a continued 
fraction tableau. Hence the conjugate of any preperiod is a 'postpcriod', an absurd 
notion. It follows that, if periodic, the two conjugate expansions are purely periodic. 

Denote by A the polynomial part of Y , and recall that Y + Y = T . It happens 
that line of the continued fraction expansion of Y + A — T is 

(38) Y + A-T = 2A-T-{Y + A-T) 

and is symmetric. In general, if the expansion of Yq has a symmetry, and if the 
continued fraction expansion is periodic, its period must have a second symmetry* . 
So if Y is exceptional in having a periodic continued fraction expansion then its 
period is of lenght 2s and has an additional symmetry of the first kind Ps = Ps+i , 
or its period is of length 2s + 1 and also has a symmetry of the second kind, 
Qs — Qs+i ■ Conversely, this is the point, if the expansion of Y has a second 
symmetry then it must be periodic as just described. 



*The case of period length 1 is an exception unless we count its one line as having two 
symmetries; alternatively unless we deem it to have period r = 2. 
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6.4. Units. It is easy to apply tlic Dirichlet box principle to prove that an order 
Q[lu] of a quadratic number field Q(w) contains nontrivial units. Indeed, by that 
principle there are infinitely many pairs of integers (j},q) so that \quj — p\ < 1/q, 
whence — (w + Lj)pq + wZUg^j < (w — ZU) + 1 . It follows, again by the box 
principle, that there is an integer I with < |/| < (oj — tU) + 1 so that the equation 
p'^ — {u! + uj)pq + iiMc^ = I has infinitely many pairs (p, q) and (p', g') of solutions 
with p = p' and q = q' (mod I) . For each such distinct pair, xl = pp' — ujujqq', 
yl = pq' — p'q + (cj + oj)qq' , yields {x — ujy){x — uJy) = 1 . 

In the function field case, we cannot apply the the box principle for a second time 
if the base field F is infinite. So the existence of a nontrivial unit x{X) —y{y)Y{X) 
is exceptional. This should not be a surprise. By the definition of the notion 'unit', 
such a unit u{X) say, has a divisor supported only at infinity. Moreover, m is a 
function of the order ¥[X, Y] , and is say of degree to, so the existence of u implies 
that the class containing the divisor at infinity is a torsion divisor on the Jacobian 
of the curve H30() . The existence of such a torsion divisor is of course exceptional. 

Suppose now that the function field ¥(X,Y) does contain a nontrivial unit u, 
say of norm —k and degree m. Then degijjY — x) = —to < — degy, so x/y is a 
convergent of Y and so some Q is ±k, say Qr = k with r odd. That is, line r of 
the continued fraction expansion of y + ^4 — T is 

line r: Yr {Y + A - T)/k = 2A/k - (Y + A-T)/k- 

here we have used the fact that (Y + Pr)/K is reduced to deduce that necessarily 

— Pj^-i^i — A — T . 

By conjugation of the (r + l)-linc tableau commencing with H38|) we see that 

line 2r: Y2r '-^ Y + A - T = 2A - T - {Y + A - T) , 

so that in any case if Y + A — T has a quasi-periodic continued fraction expansion 
then it is periodic of period twice the quasi-period. This result of Tom Berry |5] 
applies to arbitrary quadratic irrationals with polynomial trace. Other elements 
{Y + P)/Q of ¥{X, Y) , with Q dividing the norm {Y + P)(Y + P) , may be honest- 
to-goodness quasi-periodic, that is, not also periodic. 

Further, if k 7^ —1 then r must be odd. To see that, notice the identity 

(39) B[CaQ , Bai , Ca2 , Ba^ ,...]= C[Bao , Cai , Ba2 , Ca^ ,...], 

reminding one how to multiply a continued fraction expansion by some quantity; 
this cute formulation of the multiplication rule is due to Wolfgang Schmidt 
The 'twisted symmetry' occasioned by division by n , equivalent to the existence of 
a non-trivial quasi-period, is noted by Christian Friesen W. 

In summary, if the continued fraction expansion of Y is quasi-periodic it is 
periodic, and the expansion has the symmetries of the more familiar number field 
case, as well as twisted symmetries occasioned by a nontrivial n . 

One shows readily that if x/y = [A , ai , . . . , a^-i ] and x — Yy is a unit of the 
domain ¥[X, Y] then, with a^-i = kci , 0^-2 = 0,2/1^, a^-s = Ka^ , . . . , 

[ 2A - T , ai , . . . , ar-i , {2A - T)/k , a^^i , . . . , ai ] 

is the quadratic irrational Laurent series Y + A — T . Alternatively, given the 
expansion of Y + A — T , and noting that therefore deg Qr = 0, the fact that the 
said expansion of x/y yields a unit follows directly from H36|l . 
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7. Comments 

7.1. According to Gauss (Disquisitiones ArithmetictE, Art. 76) . . . veritates ex 
notionibus potius quam ex hauriri debebant^ . Nonetheless, one should not underrate 
the importance of notation; good notation can decrease the viscosity of the flow to 
truth. From the foregoing it seems clear that, given Y'^ = + 4:v{X — w) , one 
should study the continued fraction expansion of Z — ^{Y + A) , as is done in . 
Moreover, it is a mistake to be frustrated by minimal models V'^ + UV — vV = 
U^-fU + vwU. 

Specifically, we understand that - &vV = - (4/ - 1)U^ + 8v{2w - l)U 
yields Y^ ^ (A2 + 4/-1)2 + 4-8w(X-(2w-1)) by way of 2U = + Y +{Af -I) 
and {V — 8v) = XU . Now X ^ 2X + 1, Y ^ AY means that, instead, we obtain 
Y"^ = {X"^ + X + ff + Av{X - (w - 1)) . This derives from + UV ~ vV 
U'^ - fU + vwU by taking 2U = X^ + X + Y + f and V - v = XU . 

7.2. The discussion above may have some interest for its own sake, but my primary 
purpose is to test ideas for generalisation to higher genus g . An important difhculty 
when (7 > 1 is that partial quotients may be of degree greater than one without 
that entailing periodicity, whence my eccentric aside at page 1751 Happily, the 
generalisation to translating by a point (wq, Bq — ei) on the quartic model effected 
above also is a simplification in that one surely may always choose a translating 
divisor so as to avoid meeting singular steps in the continued fraction expansion. In 
that context one finds that the sequence (. . . ,2,1, 1, 1, 1, 1, 1, 2, 3, 4, 8, 17, 50, . . . ) 
satisfying the recursion ThsTh+s = Th-2Th+2 + T^ arises from adding multiples of 
the divisor at infinity on the Jacobian of the curve Y"^ = {X^ — 4X+ 1)^ + 4(Ar — 2) 
of genus 2 to the divisor [((^, 0), (^, 0)] ; ip is the golden ratio. 
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